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ABSTRACT 

Stars form out of molecular gas and supply dust grains during their last evolutionary 
stages; in turn hydrogen molecules (H2) arc produced more efficiently on dust grains. 
Therefore, dust can drastically accelerate H2 formation, leading to an enhancement 
of star formation activity. In order to examine the first formation of stars and dust 
in galaxies, we model the evolution of galaxies in the redshift range of 5 < z < 20. 
In particular, we focus on the interplay between dust formation in Type II supernova 
ejecta and H2 production on dust grains. Such effect causes an enhancement of star 
formation rate by an order of magnitude on a timescale (~ 3-5 galactic dynamical 
times) shorter than the Hubble timescale. We also find that about half of the radiative 
energy from stars is reprocessed by dust grains and is finally radiated in the far 
infrared (FIR). For example, the typical star formation rate, FIR and ultraviolet (UV) 
luminosity of a galaxy with Mvir = 10^^'^ Mq (virial mass) and Zvir = 5 (formation 
redshift) are 3 Mq yr~^, 4 x 10^ Lq, and 3 x 10^ Lq, respectively. This object will 
be detected by both ALMA and NGST. Typical star formation rates and luminosities 
(FIR, UV and metal-line luminosities) are calculated for a large set of (Afvir, ^vir)- 
Using these results and the Press-Schechter formalism, we calculate galaxy number 
counts and integrated light from high-redshift (z > 5) galaxies in sub-millimetre and 
near-infrared bands. We find that: i) ALMA can detect dust emission from severalx 10^ 
galaxies per square degree, and ii) NGST can detect the stellar emission from 10^ 
galaxies per square degree. Further observational checks of our predictions include the 
integrated flux of metal (oxygen and carbon) lines; these lines can be used to trace the 
chemical enrichment and the gas density in early galactic environments. We finally 
discuss possible color selection strategies for high-redshift galaxy searches. 

Key v^rords: dust, extinction — galaxies: evolution — galaxies: high-redshift — 
infrared: galaxies — submillimetre 



1 INTRODUCTION 



In order to understand the chemical and thermodynamical 
state of the interstellar medium (ISM) of primeval galax- 
ies, dust formation needs to be considered. Even in metal 
poor gala:x;ies, dust grains can drastically accelerate the for- 
mation rate of molecular hydrogen (H2), expected to be the 
most abundant molecule in the ISM (Hirashita, Hunt, & Fer- 
rara 2002a). Hydrogen molecules emit vibrational- rotational 
lines, thus cooling the gas. This process is particularly im- 
portant to understand the formation of stars in metal-poor 
primeval galaxies (e.g., Matsuda, Sato, & Takeda 1969; 
Omukai & Nishi 1998; Nishi & Susa 1999; Bromm, Coppi, 
& Larson 2002; Abel, Bryan, & Norman 2002; Nakamura 
& Umemura 2002; Kamaya & Silk 2002; Ripamonti et al. 
2002). The important role of dust on the enhancement of 
H2 abundance is also suggested by observations of damped 



Lya systems (DLAs; Ge, Bechtold, & Kulkarni 2001; cf. Pe- 
titjean, Srianand, & Ledoux 2000). 

The existence of dust in young galaxies is naturally ex- 
pected because Type II supernovae (SNe II) are shown to 
produce dust grains (e.g., Dwek et al. 1983; Moseley et al. 
1989; Kozasa, Hasegawa, & Nomoto 1991; Todini & Ferrara 
2001). Since the lifetime of SN II progenitors (massive stars) 
is short, SNe II are the dominant production source for of 
dust grains in young (< 1 Gyr) star- forming galaxies. The 
winds of evolved low-mass stars contribute to dust formation 
considerably in nearby galaxies (Gehrz 1989), but the cos- 
mic time is not long enough for such stars to evolve at high 
redshift {z > 5), when all galaxies should have ages smaller 
than ~ 1 Gyr. However, dust is also destroyed by SN shocks 
(McKee 1989; Jones, Tielens, & HoUenbach 1996). The de- 
tailed modelling of dust evolution in galaxies therefore re- 
quires an accurate treatment of both types of processes (for 
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recent modeling, see e.g., Edmunds 2001; Hirashita, Tajiri, 
& Kamaya 2002b). 

Here we model the evolution of dust content in primeval 
galaxies. We adopt the results of Todini & Ferrara (2001) 
for the dust formation rate in SNe H. Although further dis- 
cussion on their application of nucleation theory is neces- 
sary (e.g., Frenklach & Feigelson 1997), their results have 
been successfully applied to the interpretation not only of 
the properties of SN 1987A but also of the FIR properties 
of the young dwarf galaxy SBS 0335—052 (Hirashita et al. 
2002a). 

One of the most direct observational constraints for the 
evolution of dust content in galaxies comes from the far- 
infrared (FIR) properties of galaxies. Dust grains absorb 
stellar light and reemit it in FIR. Recent observations by the 
Submillimetre Common-User Bolometer Array (SCUBA) 
and the Infrared Space Observatory (ISO) have made it pos- 
sible to study galaxy evolution in the FIR band up to 2 < 3 
(Small et al. 1998). The detection of the cosmic infrared- 
submillimetre (sub-mm) background by the COsmic Back- 
ground Explorer (COBE) (Puget et al. 1996; Fixsen et al. 
1998) has also provided crucial information on the star for- 
mation history of galaxies in the universe (e.g., Dwek et al. 
1998). Some theoretical works have modelled the FIR evo- 
lution of galaxies up to z ~ 5 (Tan, Silk, & Balland 1999; 
Pel, Fall, & Hauser 1999; Takeuchi et al. 2001a; Xu et al. 
2001; Pearson 2001; Totani & Takeuchi 2002), and the FIR 
luminosity of galaxies per unit comoving volume seems to 
be much higher at z ~ 1 than at 2 ~ (see also Elbaz et al. 
2002). However, such a strong "evolution" beyond z — 2 has 
been excluded (Gispert, Lagache, & Puget 2000; Malkan & 
Stecker 2001; Takeuchi et al. 2001a). 

Although there is clear evidence for the existence of dust 
in galaxies at z < 5 (Armus et al. 1998; Soifer et al. 1998), 
few works focusing on the early dust formation in galaxies 
exist. Some "semi-analytic" works have included the dust 
formation in the early galaxy evolution (e.g., Devriendt & 
Guiderdoni 2000; Granato et al. 2000), but there has been 
no study treating the dust formation, the molecular forma- 
tion on grain surfaces, and the star formation history in a 
consistent manner. Therefore, in this paper, we model the 
three processes consistently so that we can obtain an obser- 
vational strategy under a consistent scenario for the early 
evolution of galaxies. 

In order to understand which physical processes gov- 
ern dust formation, observations at sub-mm wavelengths 
(300 ^m < A < 1 mm) are crucial. For high-redshift ob- 
jects, redshifted FIR radiation, i.e., sub-mm light, should 
be observed to detect the dust emission. In particular, de- 
tecting the sub-mm radiation from galaxies at z > 5 re- 
quires more sensitive and high-resolution observations (e.g., 
Takeuchi et al. 2001b). A future ground-based interferomet- 
ric facility, the Atacama Large Millimeter Array (ALMAQ), 
can be used to study such high-redshift galaxies. The de- 
tected amount of metals and stars can be used to constrain 
the galajcy evolution through a chemical evolution model 
(Tinsley 1980). Redshifted sub-mm metal emission lines can 
also be observed with ALMA. This can directly constrain 
the abundance of metals formed in the early epoch of galaxy 
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evolution. (Oh et al. 2002 have also proposed to probe high- 
redshift intergalactic medium metallicity by metal absorp- 
tion lines.) In order to detect the stellar light from the high- 
redshift universe, observations by the Next Generation Space 
Telescope (NGST^) in near infrared (NIR) will be partic- 
ularly suitable. At 2 /im, for example, we can observe the 
~ 2000 A ultraviolet (UV) light radiated from a galaxy at 
z ~ 10. Therefore, our scenario will become testable in the 
near future. This means that it is worth constructing a con- 
sistent model for the high-redshift galaxy evolution. 

For any observational facility, statistical properties of 
galaxies should be discussed to obtain a general picture of 
galaxy evolution. Two quantities are particularly important 
for statistical purposes: galaxy number counts (the number 
of galaxies as a function of observed flux) and integrated 
light (the sum of the flux from all the galaxies considered; 
Hauser & Dwek 2001 for a review). In this paper, there- 
fore, we estimate the contribution of high-redshift galaxies 
to these two quantities. 

Throughout this paper, we assume a flat cold dark mat- 
ter (CDM) cosmology with a cosmological constant. The val- 
ues of quantities are the same as those in Mo & White (2002) 
(Om = 0.3, Qa = 0.7, and Ho = WOh km s"^ Mpc"^ = 70 
km s~^ Mpc^^). The baryon density parameter is assumed 
to be f^b = 0.02/i~^. For the power spectrum of the den- 
sity fluctuation, n = 1 and erg = 0.9 are adopted. We flrst 
model the physical state of gas and the content of dust and 
metals in a galaxy during its early evolutionary stage (§ 
There we also model the luminosities of FIR, UV, and metal 
lines. The result of our model for fiducial galaxies are pre- 
sented in § ^ Based on these results, we next calculate the 
galaxy number counts and the integrated light at various 
wavelengths in § ^. We discuss the observational implica- 
tions of our results in 5 n. 



2 GALAXY EVOLUTION MODEL 

We aim at modelling the star formation rate during the early 
stage of galaxy evolution. We consider that stars form only 
in molecular environments, where cooling to a low tempera- 
ture is possible. We solve the evolution of molecular content 
in a way consistent with the star formation history and the 
dust supply from stars. Some of the past studies (e.g., Nor- 
man & Spaans 1996) have also considered similar processes, 
but we now tie star formation to the molecular content in 



determining the star formation rate (§ 2.1). 

Before z = 5, all the galaxies are younger than the cos- 
mic ag 60 of 1.2 Gyr. Although dust can play an important 
role in forming hydrogen molecules, the model for the evo- 
lution of dust content in the early epoch of the universe has 
not been developed yet. Since planned facilities will enable 
us to observe dust emission from z > 5, it is worth con- 
structing a theoretical framework for the dust evolution at 
z > 5. For observational strategies, we estimate the radiative 
energy at various wavelengths. For z > 20, the number of 
galaxies which survive stellar feedback is negligible (Ciardi 
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^ The cosmic ages for the adopted cosmology at 2 = 5, 10, and 
20 are 1.2, 0.47, and 0.18 Gyr, respectively. 
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et al. 2000) . Therefore, we concentrate ourselves on galaxies 
formed between z = 20 and 5. 



2.1 Interstellar medium evolution 

2.1.1 Gas density 

The star formation process is affected by the physical state 
of the ISM. In particular, cooling by molecular hydrogen 
plays an important role in star formation. Although the law 
governing the galactic star formation rate is not established 
at all for primeval galaxies (Kennicutt's empirical law is de- 
rived from observations of nearby galaxies; Kennicutt 1998) , 
the abundance of H2 should be ultimately a key parameter 
as stars are only seen to form in molecular complexes (see 
e.g., Wilson et al. 2000 for a recent observation of a nearby 
galaxy). Moreover, enhancement of molecular gas formation 
is shown to result in active star formation (for recent re- 
sults, see e.g., Walter et al. 2002). We define the molecular 
fraction of hydrogen, /h2 , as 



/h2 = 2nH2 / riH , 



(1) 



where tihj and wh are the number densities of hydrogen 
nuclei and hydrogen molecules, respectively (i.e., if all the 
hydrogen nuclei are in the molecular form, /hj ~ !)• The 
i mpo rtance of /hj in the star formation law is modeled in 



2.2 



We consider a pure hydrogen gas, and neglect helium in 
this paper for the following two reasons. First, helium has 
little influence on the electron abundance (i.e., the formation 
rate of H2 in the gas phase is not affected by the presence 
of He; Hutchings et al. 2002). Second, the existence of he- 
lium does not affect the temperature of gas and thus it does 
not change the reaction rates concerning H2 formation Ki- 
tayama et al. (2000). In order to examine the time evolution 
of , we solve a set of reaction equations describing the 
H2 formation and destruction. 

First, we must estimate the number density of the hy- 
drogen gas riH, because it affects both the reaction rate and 
the cooling rate. The objects that can form stars continu- 
ously against the heating from interstellar UV (lUV)Q radia- 
tion and stellar kinetic energy input have virial temperature 
typically larger than 10'* K (Ciardi et al. 2000). Since we are 
interested in the star formation activity of galaxies, we only 
consider objects whose virial temperature is larger than such 
value. If the gas temperature is larger than 10* K, cooling by 
hydrogen atomic lines is efficient. Then the gas loses pres- 
sure and collapses in the gravitational field of dark matter. 
Even in the presence of UV background radiation, objects 
with such a high virial temperature can collapse because of 
efficient cooling and H2 self-shielding (Kitayama et al. 2001). 
The collapse increases the H2 column density further, and 
as a result, external UV shielded by H2 has little effect on 
the temperature and dynamics of the gas after the collapse. 
Therefore, we neglect the external UV radiation field. 

Since the gas cooling time is much shorter than the 
Hubble timescale for the objects of interest in this paper 



■* In order to avoid confusion between external (background) and 
internal (interstellar) UV fields, we call the UV from stars within 



(Madau, Ferrara, & Rees 2001), we expect that a signifi- 
cant fraction of baryons finally collapses in the dark matter 
potential. If the halo is rotating, the gas will collapse in a 
centrifugally supported disk. The radius (rdisk) and the scale 
height (H) of the disk are determined following Ferrara, Pet- 
tini, & Shchekinov (2000) (see also Norman & Spaans 1996; 
Ciardi & Loeb 2000) . By considering the conservation of an- 
gular momentum and assuming a typical value for the spin 
parameter (A — 0.04; Barnes & Efstathiou 1987; Steinmetz 
& Bartelmann 1995), we obtain raisk — 0.18rvir, and the ra- 
dius of the dark halo, r^ii, is estimated in terms of the mass 
of the dark halo, Mvir, and the redshift of virialisation, Zvir, 
as 



— TvirAc Pc0f2M(l + Zvir)^ = Mvir , 



(2) 



where pco = Sc^ Hq /8ttG and Ac are the critical density 
of the universe at 2; = and the overdensity of an object 
formed at z^ii, respectively. The spherical collapse model 
(e.g., Peebles 1980) predicts Ac(zvir) ^ 180. We adopt the 
fitting formula given by Kitayama & Suto (1996) for Ac(2vir) 
as 



Ac(2vir) ^ 187r^[l -I- 0.4093(l/!^v - l)"-'"'^^] , 
where f2v is the density parameter at Zvir given by 

nM(l +2vir)^ 



^m{1 + 2vir)=' -f (1 - - nA)(l + 2vir)2 + S^A 



(3) 



(4) 



This formula is applicable only to a fiat universe (S^m + ^^a = 
1) with SIm < 1. 

Using Tdisk and H (the scale height of disk), nu is esti- 
mated from 



(5) 



where mn is the mass of a hydrogen atom. We have assumed 
that the mass of a galaxy is dominated by the dark halo and 
that the mass of the gas is Mvir^^b/f^M- 

An object collapsed at 2vir with a mass of Myir has a 
virial temperature, Tvir, defined as 

Tvir = ^7 , (6) 

where G is the gravitational constant, fce is the Boltzmann 
constant, and fi is the mean molecular weight. The factor 
"3" can be different depending on the radial profile of gas. 
Therefore, the definition of Tvir contains an uncertainty of 
the order unity. Since we are considering a pure hydrogen 
gas, fi = [1 + x)~ mn, where x is the ionisation degree 
(fraction of hydrogen nuclei in the ionized state) which is 
derived from equation (^ below. The circular velocity, Vc, 
is defined as 



GMvir 



we also define the circular timescale, tdr, as 
_ 27rrdiak 



(7) 



(8) 



the galaxy "lUV" 



The disk thickness relative to the radius is estimated 
by using the formalism in Ferrara et al. (2000). Using Tvir 
(eq. ^) for the gas temperature, H/r^isk is estimated to be 
~ 0.1. We fix this thickness for the disk throughout the 
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time evolution of each galaxy. Probably H could decrease 
further as gas cools, but it is difficult to include a dynamical 
evolution of gas into our one-zone model. Moreover, turbu- 
lent energy supplied by SNe 11 increases H, an effect which 
is hard to quantify in the framework of this paper. The 
resulting gas density roughly follows the scaling relation: 
riH ~ 80[(1 + 2;vir)/10]^ cm"^. This is higher than the typi- 
cal density of the local interstellar medium (wh ~ 1 cm~^), 
reflecting the high-density environment of high-redshift uni- 
verse. However, that density is lower than that calculated 
by Norman & Spaans (1996) because of our higher H/vdisk- 
Even in our case, the gas cools to reach 300 K on a short 
(~ tcir) timescale and star formation becomes possible any- 
way. The detailed treatment of H/r disk does not affect any 
of the following conclusions except for the metal-line emis- 
sion, for whi ch we will discuss the uncertainty caused by 
this factor (§ 2.5.2). The gas consumption into stars is also 
neglected. This is a valid assumption in this paper, because 
in our calculation star formation history is only traced until 
< 30% of the gas content is consumed. 



2.1.2 Chemistry and cooling 

We calculate the time evolution of ionisation degree (x), 
molecular fraction (fu^ equation |^) and gas temperature 
(T) of the hydrogen gas. The molecular fraction is particu- 
larly relevant here, because it determines the final cooling 
rate of gas when stars form (T ~ 300 K). In order to de- 
termine , we also calculate T and x. The equations are 
basically the same as those of Hutchings et al. (2002) ex- 
cept for the photo-processes (reaction of gas with photons 
from stars) and the formation of H[2 on dust grains. We sum- 
marise the reactions considered in this paper and their rate 
coefficients {Rn', n = 0, 11) in Table ^. In Table |^, we 
also list the reaction coefficient for the II2 formation on dust 
(-Rdust; see Appendix ^ for details; see also Haiman, Rees, 
& Loeb 1996 for the list of important reactions). In Table 
^ we list the photo-processes and their cross sections. The 
reaction rates of those photo-processes are expressed as r„ 
{n=12, 15). 

The time evolution of the ionizing degree is described 

as 



^ = xfoRiriH — x^i?2nH -I- ri2/o , 
at 



(9) 



where /o = 1 ~ x — /hj is the neutral fraction of hydro- 
gen. The terms on the right-hand side are the rates of colli- 
sional ionisation, recombination, and photo-ionisation. Next 
the time evolution of the molecular fraction is written as 

^ 2/oa;nH(-Roff,i + -Roff,2) + 2Rdust'DnHfo 

- fii2nii{x^RcS,3 + foRio + xRii) - Fia/Ha , (10) 



where 

Refi,! = 



R3R4 



foRi + xRs + ri4/nH 



R, 



-off, 2 



R6R7 



(11) 



(12) 



foRr + xRs + Fis/nn ' 
are the efi^ective formation rates of II2 including the efi^ect 



of destruction rate of H and Hj, respectively, and 
RgRs 



RcB,:. 



foRj + xRs + Tis/riH 



(13) 



is the destruction of due to H~ collisions. On the right- 
hand side in equation (^o|), we estimate the rates of the 
formation in the gas phase, destruction via reactions 9-11 
of Table 0, photo-dissociation, and formation on grains {V 
is the dust-to-gas mass ratio). We also list the II2 forma- 
tion rate on grains, -Rdust in Table ^ (the details are given 
in Appendix M) . The treatment of photo-ionisation and H2 
photo-dissociation rates is described in the next paragraph. 
For the initial conditions on x and /h2 a.tt = Q {z = Zvir), we 
assume the equilibrium values determined by Tvir and nn. 
The result is, however, insensitive to the choice of these val- 
ues, and as gas cools the molecular fraction always reaches 
/h2 '~ a few X 10"'^ hy t ^ 10^ yr, when the electron abun- 
dance becomes too small to produce further H2 in the gas 
phase. 

Equation (A20) of Kitayama & Ikeuchi (2000) gives F12 
(Fhi in their notation) as a function of the incident UV 
intensity and the H I column density. Although derived for 
background UV radiation, their formula is applicable to the 
lUV field as well. We assume that the photon paths are 
optically thin against ionisation of H~ and dissociation of 
HJ. We also use their equation (A21) to estimate the photo- 
ionisation heating. The effect of radiative transfer is included 
for photons ionizing H in the form of the column density 
of H I, A^'hi. We estimate the column density by Ahi — 
nH(l — X — /H2)'"disk- In order to use the formulation by 
Kitayama & Ikeuchi (2000), we describe the spectrum of 
the incident lUV radiation from stars by a power law with 
an index a: 



(14) 



where u is the frequency of photons, Io{iAii) is the inten- 
sity at the ionisation frequency of neutral hydrogen (i^hi ~ 
3.3 X 10^^ Hz). We use the same spectrum to estimate 
F„ (71 — 13, 14, 15) by extending the spectrum down to 
1.8 X 10^* Hz corresponding to the photon energy (hu) of 
0.74 eV. This corresponds to the threshold energy for the 
photo-ionisation of H~ (Table /iuv('^) could be obtained 
from a synthetic spectrum of stellar populations, but in this 
paper we simply set a = 5 for the following two reasons: 
i) The spectral shape is uncertain because of the interstel- 
lar dust extinction. In particular, little is known about ex- 
tinction curve in an extremely metal-poor environment; ii) 
Even in the case of a = 5, where the largest number of H2 
dissociating photons are produced among the four spectra 
examined in Kitayama et al. (2001), dissociation of H2 is 
negligible in the presence of dust grains. The normalisation 
of the intensity is determined from 



Luv,oexp(— Tdisk) 



r{v) dv , 



(15) 



where we define !/min as the minimum frequency where OB 
stars dominate the radiative energy of star-forming galaxies 
(~ 10^^ Hz), I/uv,o is the UV luminosity estimated in equa- 
tion (p^), and Tdisk is the typical dust optical depth in the 
disk. This typical optical depth can be simply estimated by 
multiplying the typical column density of dust in the direc- 
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Table 1. Reaction rates needed to calculate the abundance of H2. For the references, see 1) Omukai (2000); 2) Galli & Palla (1998); 3) 
this paper. The unit of the gas temperature T is K unless otherwise stated 



No. 



Reaction 



Rate [cm^ 



Ref. 



1 


H4 


- e — 


.H+ + 


2e- 


exp[-32.71 4- 13.54 ln(r(eV)) - 5.379(ln(T(eV)))2 4- 1.563(ln(T(cV)))3 




1 












- 0.2877(ln(T(eV)))4 4-3.483 x 10-2(ln(r(eV)))5 - 2.632 x 10-3(ln(T(eV)))'' 














4- 1.120 X 10-*(ln(T(eV)))^ - 2.039 x 10-«(ln(r(eV)))8] 






2 


H+ 


4-e- - 


H 4- 


7 


exp[— 28.D1 — 0.7241(ln(_i (eV))) — 2.0zd x 10 (ln(J (eV))) 
- 2.381 X 10-3(ln(T(eV)))3 - 3.213 x 10~"'(ln(r(eV)))-' - 1.422 x 10" 
4- 4.989 X 10-'3(ln(T(eV)))'3 4- 5.756 x 10-'^(ln(r(eV)))7 - 1.857 x 10" 
-3.071 X 10-''(ln(T(eV)))9] 


5(ln(r(eV)))5 
8(ln(r(eV)))8 


1 


3 


H4 


- e~ — 


^ H- 4- 


7 


1.4 X io-i8rO-928 exp(-T/1.62 X 10^) 




1 


4 


H- 


+ H - 


-^H2 4 


e— 


4.0 X lo-^T-"-!'' (T > 300); 1.5 x 10^^ (T < 300) 




1 


5 


H- 


+ H+ 


— > 2H 




5.7 X lO-'^T-i/^ 4- 6.3 x 10-* - 9.2 x lO-^T^/^ + 4.4 x lO-^^T 




1 


6 


H4 


-H+ - 


-+H+4 


-7 


dex[-19.38 - 1.523 login ^ + l-118(logio T)'^ ' 0.1269(logi(, T)3] 




1 


7 




+ H - 


-^H2 4- 


H+ 


6.4 x 10""' 




1 


8 




4-e- - 


2H 




2.0 X lO-^T-1/2 




1 


9 


H2 + H+ - 




+ H 


3.0 X 10-i°exp(-21050/r) (T < 10'') 




2 












1.5 X 10-i''exp(-14000/r) (T > 10'') 






10 


H2 


+ H — 


^ 3H 




H L 

A;l = 1.12 X 10-10 exp(-7.035 x lO^/T) 

ku = 6.5 X lO-^T-1/2 exp(-5.2 x 10*/T)[1 - exp(-6000/T)] 
a = 4.0 - 0.416 logio(T/10'*) - 0.327(logio(T'/10''))2 




1 


11 


H2 


4-e- - 


2H 4 


- e— 


4.4 X 10-i°T°-35exp(-1.02 x 10^ /T) 




1 


dust 


H4 


' H + grain — 


* H2 4- grain 


2.8 x 10-i'5(T/100 K)i/2 if T < 300 




3 












if T > 300 see Appendix ^ 







Table 2. Cross sections for photo-ionisation and photo-dissociation processes. The unit of u (frequency of light) is Hz. For the references, 
see 1) Kitayama & Ikeuchi (2000); 2) Abel et al. (1997); 3) Tegmark et al. (1997) 



No. 


Reaction 






cross section 
(cm^) 




V range 
(Hz) 


Ref. 


12 


H + 7 . 


H+4 


- e 


6.30 x 10-i**(;//3.3 x lO^^)-^ " 




!^ > 3.3 X lO^s 


1 










effect of optical depth also included 


in the text 






13 


H2 4-7 — 


-H* 


— > 2H 


see equation (jl^ 






2 


14 


H- -H7 - 


H4 


- c^ 


3.486 X 10-1^ (a' - l)3/2/x='-" {x = 


vl\.% X 10") 


!^ > 1.8 X 10^'' 


3 


15 


H++7- 


H4 


-H+ 


7.401 X 10-1** dox(-a;2 - 0.0302x3 _ 


- 0.0158x4) 


1/ > 6.4 X lO^" 


3 










(x = 2.7621n(i//2.7 X lO^^) 









tion of disk plane, ndust^dust, by the absorption cross section 
of dust against UV light, tto^Quv, as 



Tdisk = 7ra ndustJ'disk • 



(16) 

where a is the grain radius (spherical grains are assumed). 
The meth od for estimating Luv.o and ndust will be described 
later in 5 2.5.1. 



The H2 photo-dissociation cross section is estimated 
from the rate given by Abel et al. (1997). However, if 
the column density of H2 becomes larger than 10^'' cm"'^, 
self-shielding effects become important (Draine & Bertoldi 
1996). Therefore, we use the following expression for the H2 
dissociation rate: 



ri3 



(47r)l.l X 10*7iuv(3.1 x 10^^ Hz) 
'Wh, /Ha ''disk' 



/ WH/Ha^diskA ' 
V 10" cm-2 ) 



(17) 



where 7iuv is in cgs units. 

For the H2 formation on grains, recent experimental re- 



sults have indicated that S ~ (the sticking efficiency of 
hydrogen atoms; Appendix ^ for Tdust > 20 K (Katz et 
al. 1999). Such a low-temperature threshold for the H2 for- 
mation suggests the following scenario. Because of the ther- 
mal coupling between the dust and the cosmic microwave 
background (CMB) of temperature Tcmb = 2.7(1 + z) K, 
molecular formation on grains might have been strongly sup- 
pressed when Tcmb > 20 K or ^ > 7. As a result, galaxies 
start to form stars actively when z ^ T .\n this case, 2^7 
is the typical redshift for the onset of active star formation. 
Since such a typical redshift, if it exists, could be in prin- 
ciple detected as an enhancement of galaxy number counts 
at a certain flux level, it is worth considering such a typical 
"formation epoch" in this paper. We reconsider this point 
in detail in § W. 

SNe H affect the abundance of H2 as they create regions 
filled with hot gas (T ^ 10^ K). In such an environment H2 
is destroyed, but it reforms after the gas cools; this pre- 
dominantly occurs in cooled shells (Shapiro & Kang 1987; 
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Ferrara 1998; we note that our /hj is roughly two times that 
of Shapiro & Kang) and becomes a few x 10""^. We can- 
not consider these effects in this paper, because our model, 
which only treats averaged quantities over the whole galaxy, 
cannot describe local bubbles and shells produced by SNe 
II. We expect that the formation of H2 on dust surfaces and 
in the cold shells keeps /h2 from decreasing significantly. 

In order to calculate the temperature evolution, cooling 
and heating should be included in our model. We adopt the 
cooling functions summarized in § 2.3 of Hutchings et al. 
(2002), i.e., cooling by molecular hydrogen, and coUisional 
excitation and ionisation of atomic hydrogen. Again cooling 
by helium is neglected because Hutchings et al. (2002) have 
shown that the temperature evolution is little affected by the 
helium cooling. For the heating by stellar lUV radiation, we 
adopt equation (A21) of Kitayama & Ikeuchi (2000). The 
initial value for Tgas is assumed to be Tvir. 



2.2 Star formation law 

The adopted star formation law is central to this paper, 
because we propose a new "paradigm" based on the fact 
that the stars are formed during the final cooling by molec- 
ular hydrogen. If the molecular gas is abundant, stars are 
formed as a result of a dynamical collapse of gas. Since a 
representative timescale of the dynamics of the gas disk is 
the circular timescale, tcir, we expect that the star formation 
rate tp is roughly Mgas/tdv in a fully molecular gas. Most of 
the "semi-analytic" recipes of galaxy evolution (e.g., Kauff- 
mann & Chariot 1998; Somerville & Primack 1999; Cole et 
al. 2000; Granato et al. 2000; Nagashima et al. 2001; cf. 
White & Frenk 1991) assume this kind of law (the depen- 
dence on the circular velocity is nonlinear though). However, 
since the final coolant is molecular hydrogen, we should in- 
clude the hydrogen content in the formulation to obtain a 
more physical star formation law. Therefore, we assume the 
following form for the star formation rate: 



V'(i) =/H2(t) A/gas Acir, 



(18) 



where Mgas and tcir are taken to be constant in time. 

Although the following results are critically dependent 
on the assumed form of /hj, the experimental evidence 
shows that star formation depends positively on molecu- 
lar abundance (e.g., Rana & Wilkinson 1986; Wilson et al. 
2000; Walter et al. 2002; but see e.g., Buat^eharveng, & 
Donas 1989; Tosi & Diaz 1990). Equation (Q is given "as 
a first approximation" in this paper to include this experi- 
mentally supported law. Although multi-phase behavior of 
ISM can be also important (McKee & Ostriker 1977; Ikeuchi 
1988; Norman & Spaans 1996), this equation is expected to 
approximate the star formation rate of the whole galaxy in 
the framework of our one-zone treatment. 

We are interested in objects whose virial temperature is 
typically larger than lO" K (§ O). Madau et al. (2001) have 



argued that such objects cool on a timescale much shorter 
than the dynamical timescale and consequently experience 
an initial strong episode of star formation. Our approach 
in this paper is conservative in the sense that active star 
formation does not occur until a significant amount of hy- 
drogen molecules is produced. However, even in our "con- 
servative" treatment, a burst of star formation occurs as a 
result of molecular formation on dust grains (§ Therefore, 



a scenario similar to, albeit physically different from, those 
of Madau et al. (2001), Ciardi et al. (2000), and Barkana 
(2002) emerges from the present study. 



2.3 Evolution of dust content 

Because of short cosmic timescale between z = 5 and 20 
(< 1 Gyr), the contribution of Type la SN and winds from 
late-type stars to dust formation is assumed to be negligi- 
ble. In this case, SNe II are the dominant sources for dust 
formation. The rate of SNe II as a function of time, 7(t), is 
given by 



(m) dm , 



(19) 



J 8 Mq 

where i/)(f) is the star formation rate (SFR) at t (for t < 
0, V(^) ~ 0), 0(m) is the initial mass function (IMF; the 
definition of the IMF is the same as that in Tinsley 1980), 
Tm is the lifetime of a star whose mass is m, and we assumed 
that stars with m > 8 Mq produce SNe II. In this paper, 
we assume a Salpeter IMF {(l){in) oc m~^ "^^) with the stellar 
mass range of 0.1-60 Mq. It has been suggested that the 
IMF is much more weighted to massive stars, i.e., top-heavy, 
in primeval galaxies (e.g., Bromm et al. 2002). In a top- 
heavy environment, the production of dust by massive stars 
is enhanced. As a result, dust amount, and thus molecular 
amount, in galaxies would be larger than that predicted in 
this paper. Here, we "conservatively" assume the Salpeter 
IMF. 

Dust destruction by SNe II can be important. The de- 
struction timescale tsn is estimated to be (McKee 1989; 
Lisenfeld & Ferrara 1998) 



M^ 



TSN = 



7eMs(100 km s" 



(20) 



where Ms(100 km s"^) = 6.8 x 10^ PIq (Lisenfeld & Fer- 
rara 1998) is the mass accelerated to 100 km s~^ by a SN 
blast, 7 is the SN II rate, e ~ 0.1 (McKee 1989) is the efli- 
ciency of dust destruction in a medium shocked by a SN II. 
Since we are interested in the first star formation activity, 
we assume the relation between stellar mass and lifetime of 
zero-metallicity stars in Table 6 of Schaerer (2002) (the case 
without mass loss is applied). 

Then the rate of increase of Md is written as 

—r- = md7 , (21) 

at TSN 

where md is the typical dust mass produced in a SN II. To- 
dini & Ferrara (2001) showed that rud varies with progeni- 
tor mass and metallicity. There is also some uncertainty in 
the explosion energy of a SN II. The Salpeter IMF-weighted 
mean of dust mass produced per SN II for the 1) Z — 0, 
Case A,2)Z = 0, Case B, 3) Z = W'^Zq, Case A, and 4) 
Z = IO'^Zq, Case B are 1) 0.22 Mq, 2) 0.46 Mq, 3) 0.45 
Mq, and 4) 0.63 Mq, respectively {Z is the metallicity, and 
Cases A and B correspond to low and high explosion en- 
ergy|^, respectively). We adopt the average of the four cases, 
i.e., md — 0.4 Mq but we should remember that md can 



^ The kinetic energies given to the ejecta are 1.2 X 10 ergs 
and ~ 2 X 10^^ ergs for Case A and Case B, respectively. 
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have a range (from 0.22 to 0.63 Mq). Since the dust de- 
struction is negUgible as shown later, our final value of the 
dust mass is approximately proportional to the adopted one 
for ma. 

2.4 Evolution of metal content 

The evolution of metal content can be predicted once the 
star formation history and metal yield per SN II are fixed 
(e.g., Tinsley 1980). Suginohara, Suginohara & Spergel 
(1999) have proposed that the fiuctuations in space and 
wavelength of (sub-) millimetre background radiation made 
of high-redshift metal lines can be used as an indicator of 
structure formation at high redshift. Here we obtain the 
metal mass injected in the gas phase by subtracting the dust 
mass from the metal mass. The evolution of the mass of a 
heavy element (species i) in the gas phase of a galaxy, Mi, 
is thus calculated by 

ir^'"'^-^^' (22) 

where rrii is the averaged mass of element i formed per SN 
II and Mdust.i is the mass of element i in the dust phase. 
According to Todini & Ferrara's calculation, dust contains 
15% of oxygen and 36% of carbon. Therefore, we assume 
that Mdust.o = O.lSMdust and that Mdust.c = 0.36Mdust 
(we take the same average as that in § 2.3). In Table ^ 
we list rrii for carbon and oxygen. When calculating rrii, we 
adopted the results by Woosley & Weaver (1995) and took 
the same mean as that in § 2.5 . Suginohara et al. (1999) also 
considered nitrogen, but we have not included this species 
as the production of nitrogen in a SN II is one or two orders 
of magnitude less than that of oxygen or carbon. 



2.5 Radiative properties 

2.5.1 UV and FIR 

Probably the most direct way to reveal high-redshift dust 
is to observe FIR emission. We now derive the evolution of 
galactic FIR luminosity and dust temperature. Because of 
the large cross section of dust against lUV light and the 
intense lUV radiation field in a star-forming galaxy, we can 
assume that the FIR luminosity is equal to the absorbed 
energy of lUV light. 

First, we should estimate the fraction of the lUV radia- 
tion absorbed by dust. For the convenience of the following 
calculation, we define the following typical optical depth in 
the vertical direction of the disk, tq, as 



To = vra Quv".dust-ff . 



(23) 



In this paper we assume single value for a, because there 
is a typical size of dust produced by SNe II as shown by 
Todini & Ferrara (2001). We note that this optical depth is 
different from Tdust in equation (p^, where we needed an 
optical depth in the disk direction. Here, the optical depth 
in the vertical direction is useful as we see in the following. 

Since tq is independent of H , the optical depth of the 
dust does not depend on the treatment of rdisk/H in § 2.1.1 . 



The dust density is related to the mean dust number density 
as 

^a^5ndu.tvrr^isk2i^ = Mdu.t , (24) 



where 5 is the grain material density. By solving equation 
( p^ for Tidust and substituting it into equation (^3|), we ob- 
tain 



To = 



3 QvvMi 



(25) 



If the angle between the direction of a lUV light propagation 
and the vertical direction of the disk is 6, the intensity of 
the light becomes roughly exp(— ro/ cos S). Therefore, the 
luminosity of lUV light escaping from the galactic disk, Luv 
is estimated to be 



■,o ^exp 



To 



LvY,0 E2{to) , 



(26) 



where {■)g indicates the mean over the solid angle and Luv.o 
is the intrinsic UV luminosity of the galaxy. The exponential 
integral Enix) (n — 0, 1, 2, a; > 0) is defined as 



E„{x) 



exp(-y) 

y 



dy . 



(27) 



Luv.o is assumed to be equal to the total luminosity of OB 
stars whose mass is larger than 3 Mq (Cox 2000): 



-At) 



dm 



dt' L{m)4>{m)ip{t~t'), (28) 



where L{m) is the stellar luminosity as a function of stellar 
mass (m). For L(m), we adopt the model of zero-metallicity 
stars without mass loss in Schaerer (2002). We fix this rela- 
tion as a first approximation in this paper. 

We adopt a ~ 0.03 /im (Todini & Ferrara 2001), Quv ^ 
1, and 5 ~ 2 g cm~^ (Draine & Lee 1984). We only consider a 
single value of a, because there is a well defined sharp peak 
in the dust size distribution by Todini & Ferrara (2001). 
Then we obtain the energy absorbed by dust. We assume 
that all the absorbed energy is reemitted in the FIR. Thus, 
the FIR luminosity Lfir becomes 



(29) 



The dust temperature Tdust is determined from the 
equation derived by Hirashita et al. (2002a) based on the 
dust emissivity given by Draine & Lee (1984) for a — 
0.03 ^im: 



Td = 20 



2.5 X lO^Md/Mo 



1/6 



K. 



(30) 



2.5.2 Metal lines 

We next calculate the luminosities of the metal lines. Our 
model calculation indicates that almost all the gas is in neu- 
tral form (x <C 1 and /hj ■C 1); hence we consider lines 
typical for neutral regions. The most interesting lines are 
listed in Table ^. For spontaneous transition from an upper 
level « to a lower level /, the total luminosity of the line with 
frequency i/unc (line = C609, C370, 063, 0146 for C I 609 
/im, C I 370 /im, O l 63.2 ^m, and O I 146 /^m, respectively) 
becomes 



(31) 



where hp is the Planck constant, Mu is the total number of 
the atoms in the upper level in the whole galaxy, and j4„i 
is the Einstein coefficient of the spontaneous emission. We 
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assume that all the carbon and the oxygen are neutral; we 
will also consider the C ii 158 /xm line, which is also emitted 
from the neutral medium (Tielens & HoUenbach 1985; Liszt 
2002). If the timescale of the spontaneous emission is longer 
than that of the coUisional excitation, the fraction within the 
state is proportional to 2g+l (i.e., 1/9, 3/9, and 5/9 for ^Po, 
^Pi, and ^P2, respectively). However, if the density is so low 
that the coUisional excitation does not occur so frequently 
as the spontaneous emission, the population in the upper 
level is reduced. We include this effect of low density by 
multiplying the factor described in Appendix ^. 

The line intensity is estimated for the case where the 
population ratio is determined by 2(7 + 1. Since this gives 
us the maximum line intensity, we denote the intensity as 
L™n^; more realistically, the line intensity is closer to !FL]^^^. 
We estimate J\fu from the metal mass in equation (|2^) and 
the number fraction in the upper level. The actual intensity 
of metal line is reduced by a factor JT. However, a precise 
determination of J- requires a detailed tre atme nt of the evo- 
lution of gas density. As discussed in § |2.1.l| , our present 
framework can only provide a first order estimate of such 
quantity; therefore, only V^^c is calculated in this paper. 
Nevertheless we propose a simple method to estimate J- in 
Appendix ^ When we calculate the C ll 158 /xm line inten- 
sity, we assume that all the carbon atoms are in the form of 
C+. 



3 RESULTS FOR FIDUCIAL GALAXIES 
3.1 A typical primeval galaxy 

In the following we show the evolution of some character- 
istic quantities predicted by our model for a galaxy with 
(Mvir, 2vir) = (10^ Mq , 10), Corresponding to a 2.5-0- den- 
sity fiuctuation of the cosmological density field. The quali- 
tative behavior of the quantities are similar for other objects 
of interest to the present study. In the next subsection, we 
explore various (Mvir, ^vir) values. 

Using the equations above, we simultaneously calculate 
the time evolution of /h2, ■0, Mdust, and various kinds of 
luminosities in a self-consistent manner. In Fig. |l|, we show 
the evolution of the molecular fraction (solid line of Fig. 
|l^), and the star formation rate ip (solid line of Fig. |^b) for 
an object Zvir ~ 10 and Mvir = 10^ Mq. We show the results 
up to lOfcir (— 7.8 X 10* yr), which is comparable to the 
time between z = 10 and z = 5. At t = lOtcir, ~ 30% of the 
gas is converted into stars. Therefore, the calculation further 
than this time is not consistent with our assumptions, which 
neglect the gas conversion into stars. The gas temperature 
rapidly drops and reaches < 300 K within lO'^ yr. 

For other sets of (Mvir, Zvir), an enhancement of the 
star formation at 3-5tcir is commonly seen. This results from 
the H2 formation on dust grains. The dotted lines of Figs. 

and |l|b indicate the case in which we neglect the dust 
(i.e., Md — 0) at all times. In this case, after stars form, 
/hj decreases because of photo-dissociation. Comparing the 
solid and dotted lines in Figs. ^ and b, we conclude that 
the existence of dust is essential in causing a strong star 
formation activity in primeval galaxies. 

In Figure |^ we show the evolution of dust, oxygen and 
carbon (solid, dotted, and dashed lines, respectively). The 
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Figure 1. Time evolution of (a) tlie molecular fraction and 
(b) the star formation rate for a galaxy with Zvir = 10 and 
Mvir = 10® Mq corresponding to the 2.5-0" density fluctuation. 
The evolution in the first 10 circular times is shown. The solid 
line is for the case where we consider dust (and therefore the H2 
formation on grains), and the dotted line is for the case of no dust 
production (A/duat = 0)- The vertical dashed line marks the 4 cir- 
cular timescales, where the typical quantities are defined in the 
text and in Table |3| In the top axis, we show the corresponding 
redshifts. 



mass fraction relative to the total gas mass is also shown on 
the right axis. Note that the solar abundances for oxygen 
and carbon are 1.0 X 10"^ and 3.4 x 10~^ in mass ratio, 
respectively (Anders & Grevesse 1989; Cox 2000). We see 
rapid rises of those masses around t ~ itdi (marked with the 
vertical dashed line). We also show the case in which dust 
destruction by SNe II is neglected (dot-dashed line): dust 
destruction has little infiuence on the dust amount during 
the early (< several tdv) evolution. 

In Figure ^ we show the evolution of three luminosi- 
ties; Luv, ipiR, and I/oi46 (solid, dotted, and dot-dashed 
lines, respectively). Whereas Lfir and ioi46 grow rapidly 
around t ~ 4tcir, L\jy does not increase as significantly as 
the other two luminosities because of dust extinction. The 
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Table 3. Typical quantities for various sets of {Afvir,2vir) 
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FIR luminosity becomes comparable to the UV luminosity 
soon after the active star formation phase starts. Therefore, 
even in the early stages of star formation, dust absorbs a 
significant fraction of stellar light and reprocesses it to the 
FIR range. This is because of the dense and compact nature 
of high-redshift galaxies (rvir oc (1 + z)~^ with a fixed Afvir). 
Therefore, even at early cosmic epochs (z > 5) simultaneous 
observations of FIR emission from dust and of the UV/NIR 
stellar light is crucial to trace the total radiative energy from 
galaxies. 

3.2 Dependence on (Mvir, Zvir) 

Table ^ clarifies the dependence of various quantities on 
Mvir and Zvir- We define the typical star formation rate for 
each object, i/>(Mvir, Zvir), as i/>(4tcir). This is because after 



4 circular times objects have consumed ~ 5% of the gas. 
Since it is empirically known that the star formation effi- 
ciency during an episode of star formation activity is < 10% 
(e.g., Inoue et al. 2000; Barkana 2002), the definition yields 
a good working value. Moreover, the value at t = 4tcir gives 
a good average for the star formation rates over the Hub- 
ble timescale at Zvir within a factor of ~ 2. In Table ^ we 
show the typical star formation rate as a function of Zvir and 
Mvir. In Table ^ we also show other quantities, all of which 
are estimated at t — Atdv (i.e., the quantities with the upper 
bars indicates the typical quantities estimated at t = 4tcir). 

All quantities except for Luv are monotonically increas- 
ing functions of both Mvir and Zvir. An object with larger 
Mvir tends to contain a larger gas mass (Mgas). An object 
with larger Zvir tends to have a short icir (oc (1 -I- Zvir)~^'^) 
because of a small rdisk oc (1 -f Zvir)~^ and a large tic oc 
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Figure 2. Time evolution of the masses of dust, oxygen (in gas 
phase), and carbon (in gas phase) in the same galaxy as Fig. |l| 
(solid, dotted, and dashed lines, respectively). The dot-dashed 
line represents the case in which dust destruction is neglected. 
The vertical dashed line shows the 4 circular timescalcs. 
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Figure 3. Time evolution of the ultraviolet luminosity (solid 
line), the far-infrared luminosity (dotted line), and the maximum 
line luminosity of O I 146 fim (igiJe) (dot-dashed line) for the 
same galaxy as in Fig. ^ The vertical dashed line shows the 4 
circular timescalcs. 



(1 -|- Zvii)^ "'. Therefore, the star formation law (eq. [l^ ) in- 
dicates that ?/j becomes an increasing function of both Afvir 
and 2vir- The larger the star formation rate is, the larger 
the metal production, the dust production, and the intrinsic 
OB star luminosity. As a result almost all quantities be- 
come an increasing function of both Mvir and Zvir- On the 
contrary, Luv is not monotonic. For objects with z^h > 15 
and Mvir > lO" Mq, Luv decreases as Mvir and Zvir in- 
crease. This is because of the compactness of high-redshift 
objects. Dust accumulation in a compact region makes the 
optical depth of dust very large; hence UV radiation is ef- 



ficiently absorbed by dust. However, objects with Zvii > 15 
and Mvir ^ 10^" Mq are extremely rare because they corre- 
spond to density fluctuations > 4 cr. 

The temperature of dust is generally higher than that 
observed in normal spiral galaxies (~ 20 K), because the 
compactness (meaning high density) of the high-redshift 
galaxies leads to a strong stellar radiation field. For a nearby 
dwarf galaxy, SBS 0335—052, the compactness of the star- 
forming region is also suggested to result in a high dust 
temperature (Dale et al. 2001; Hirashita et al. 2002a). 



4 STATISTICAL PREDICTIONS 

Future observations will provide large samples of galaxies at 
2 > 5. Therefore, predictions of the statistical properties of 
galaxies in such a redshift range are very valuable. Using the 
result above, we focus on cosmic star formation rate (Tins- 
ley & Danly 1980; Madau et al. 1996), number counts of 
galaxies, integrated light (Partridge & Peebles 1967), where 
"integrated" means the sum of the contribution from all the 
galaxies of interest. 



4.1 Cosmic Star Formation History 

The star formation rate per unit comoving volume at the 
redshift z, ^'(z), is written as 



*(z) 



dMtpiM, z 



dn(M, z) 
dM 



(32) 



where n(M, z) is the comoving number density of halos with 
masses larger than M at z, and tp{M, z) is the star formation 
rate of an object with mass M at z. 

In equation (^), Mmin is set equal to the maximum 
mass of a galaxy whose gas is blown away by SNe II in order 
to exclude objects which do not form stars continuously. 
From a fitting to the result of Ciardi et al. (2000) , we adopt 



5 X 10 
(1 + ^) 



T7. Mq. 



(33) 



Ciardi et al. presented three cases (A, B, and C) for various 
cooled gas fraction, star formation efficiency, and escaped 
photon fraction^ all of which are inherent to their model. 
In equation (pSl), we adopt their fiducial case A. Even for 
the other cases, the following results are not changed at all, 
because galaxies whose mass is as low as 10^ Mq do not 
contribute to the number counts at the flux level of interest. 
Ciardi et al.'s case A indicates that the 1.2% {= fhf*) of the 
baryons are converted into stars in a dynamical time. This 
gives a star formation rate comparable to the one that we 
find at t = 4tcir. 

For the maximum mass we adopt Afmax ~ W^"^ Mq. As 
long as Mmax ^ Mq , the results are not changed at all, 
because the number of objects with Mvir > 10^'' Mq is neg- 
ligible in all the considered redshift range (5 < z < 20). To 
derive dn{M, z)/dM, we use the Press-Schechter formalism 
(Press & Schechter 1974). We assume that the star forma- 
tion rate is a function of the virial mass and the redshift as 
listed in Table 0; i.e., ?/'(M, z) = i/'(Mvir = M , Zvir = z). 

In Figure M we show ^'(z) (solid line). The consid- 
ered range of redshift is between Zrnin ~ 5 and Zmax ~ 
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Figure 4. Star formation rate per unit comoving volume (solid 
line). The dotted line indicates the star formation history with a 
"burst" (bSFH) at z = 6.5. The two data points around z ^ A are 
from Madau et al. (1998) (lower) and Steidel et al. (1999) (up- 
per). These points have been corrected for the dust extinction by 
Steidel et al. (1999) and for the cosmological parameters chosen 
in this paper. 



20, because we are interested in the first cosmic Gyr (~ 
Hubble timescale at z — 5) and the number of galaxies with 
2max > 20 is negligible. The comoving stellar mass den- 
sity formed by z; = 10, 7, and 5 are 2.2 x lO'^ Mq Mpc"^, 
2.7 X 10'' Mq Mpc"^, and 1.3 x lO'^ Mq Mpc"^, respectively. 
In units of the critical density, these numbers correspond to 
= 1.6 x 10"'^, 2.0 X 10'^, and 9.6 x 10"^, respectively. 
As stated in '• 



2.1.2 , the suppression of H2 formation on 



grains for Tdust > 20 K as shown experimentally by Katz 
et al. (1999) will introduce a characteristic star formation 
epoch (zburst). Their results indicate that H2 formation on 
grains is allowed after the CMB temperature, Tcmb, drops 
to < 20 K (z < 7). Since the H2 formation on dust is crucial 
to activate star formation in a galaxy, galaxies start to form 
stars actively at Zburst ~ 7 in this scenario. Motivated by 
this, we have also investigated another star formation his- 
tory activated at 2: ~ Zburst- We call such a star formation 
history "bursting star formation history" (bSFH). On the 
other hand, we call the star formation history without ^burst 
"continuous SFH" (cSFH; the solid line of Fig. |). Unless 
otherwise stated, we show the results for the cSFH. 

Galaxies with an active star formation are found up to 
z ~ 6.5 (e.g., Hu et al. 2002). Moreover, Shanks et al. (2001) 
have shown that the space density of bright galaxies at z ~ 6 
is comparable to the local space density (see also Ouchi et 
al. 2002) , which suggests that the formation epoch of galax- 
ies lies at z > 6.5. In this paper, we adopt the lower limit 
ZhuTst = 6.5. As shown later on, even for this low Zburst value, 
ALMA and NGST cannot clearly identify the existence of 
Zburst- However, it is worth examining the typical flux level 
where the effect of Zburst on number counts is clearly seen. 

Even in the presence of the typical burst redshift, Zburst, 
the number of the virialised galaxies does not change. The 
only difference is that the star formation is largely sup- 
pressed for z > Zburst in the bSFH. Moreover, the star for- 



mation rate is typically enhanced by a factor of 5 at Atdi 
due to the enhancement of molecular formation on accumu- 
lated dust. Therefore, for ^'(z) in the bSFH, we adopt a 
value equal to 1/5 of the one calculated in equation (^2|) if 
z > Zburst- For z < Zburst, wc adopt the same value both for 
cSFH and for bSFH. The bSFH with Zburst = 6.5 is shown 
by the dashed line in Fig. ^. 

4.2 Galaxy number counts 

We estimate the number of galaxies with observed flux 
greater than at an observed frequency f by 

dz / dM 

dViz)' 



dn(M, z) 
dM 



dz 



(34) 



where Mu^{f^, z) corresponds to the mass of the galaxy 
with speciflc luminosity Li,, and 



(l + z)L,(i 



(1+^) 



(35) 



The luminosity distance, di^ is given by (Carroll, Press, & 
Turner 1992) 



dL =c(l + 2) / dz {1 + z') 



dt 



dz' 



(36) 



where c is the light speed, and 
dt 



dz 



Ho{l 



X ^(1 + nM2)(l + z)2 - nAz{2 + Z) . (37) 

The volume element per unit redshift, dV/dz, is written as 



dV{z) 
dz 



4ncdl 
(1+^ 



(38) 



We assume that depends on (M, z) and that mono- 
tonically increases as M for any z and any z^. The mono- 
tonicity is guaranteed in the sub-mm number counts, but 
not in the NIR counts. A NIR observation can detect the 
UV light in the restframe of galaxies, and the UV luminos- 
ity is not a monotonically increasing function of Afvir for 
Zvii > 15 (Table Therefore, in calculating the NIR num- 
ber counts, we set ^max = 15. If Mumifu, z) > Mmax, we set 

Mim(/i^) equal to Mmax- 

Finally, we must flx as a function of M and z. We 
approximately identify (M, 2) with (Mvir, Zvir), and use the 
quantities listed in Table ^ as a function of (Mvir, 2vir). The 
details of the calculation will be described separately in the 
following two subsubsections for sub-mm and NIR, respec- 
tively. 



4-2.1 Submillimetric bands 

We calculate the number counts due to continuum emission 
from dust. The peak of dust emission falls in the FIR range 
(~ 50-150 /im) in the restframe of a galaxy. The peak of 
high-redshift (2 > 5) galaxies is shifted to the sub-mm range 
in the observational restframe. For example, for a galaxy at 
2 ~ 7, observations around ~ 800 fim should be performed. 
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Figure 5. Theoretical predictions of the number counts of the high-redshift (5 < ^ < 20) galaxies in the ALMA bands (solid lines): (a) 
450 /im, (b) 850 ^m, and (c) 1.3 mm. The number counts of galaxies whose redshift range is 5 < z < 7 are also shown (dotted lines). 
The dashed lines show the result for burst mode star formation (bSFH), in which galaxies start to form stars actively at z = 6.5. The 
vertical dot-dashed lines show the ALMA detection limits. 



Here we consider three ALMA bands centered at 450, 850, 
and 1300 /im. 

The spectral energy distribution of the FIR emission 
from a galaxy is assumed to be a modified blackbody with 
a temperature Tdust (listed in Table |^ , 

= C//^B4fdu,t), (39) 

where the coefficient C is determined so that the integra- 
tion of Lv for all the wavelength range becomes equal to 

i^FIR(Mvir = M, 2vir = z) . 

From the above we can relate M and for each z 
using the relation between L„ and (Zfir, Tdust), and then 



equation ( p4| ) to obtain N{fv, v). In Figs. |5p,-c, we show the 
number counts in three of the ALMA observational bands 
of 450, 850, and 1300 ^m for galaxies with 5 < 2 < 20 
(solid lines) and with 5 < 2 < 7 (dotted lines). The vertical 
lines show the detection limits, for which we adopt the same 
values as Takeuchi et al. (2001b) for the 5 a limits of 8 hr 
integration (220, 16, and 4.6 /iJy, respectively). We see that 
5, 1.5 X 10^, and 2.9 x 10^ galaxies per square degree can be 
detected in the 450, 850, and 1300 ^m bands, respectively. 
The high angular resolution of ALMA enables us to detect 
those galaxies without confusion. A bright source can make 
the detection of faint sources difficult because of a limited 
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dynamic range of the detector. The probabihty that sources 
which are 10^ (this number comes from the dynamic range of 
the detector of ALMA) times larger than the 5 a detection 
limit of ALMA exist in the field of view is negligible (< 
3 X 10"''; Takeuchi et al. 2001b). From Figure |, it is also 
concluded that ALMA is sensitive to galaxies with z < 7. 



4-2.2 Near infrared bands 

In order to test our theoretical prediction further, it is neces- 
sary to detect the stellar light, because dust production his- 
tory is deeply related to star formation history. Since young 
galaxies are characterized by a strong UV stellar continuum 
produced by OB stars, we examine the detectability of the 
UV light, which is redshifted to the NIR range in our obser- 
vational restframe. For example, an observation at 1.8 j-im 
can detect the stellar light at 3000 A in the restframe of a 
galax;y at 2 = 5, or that at 2000 A at 2 = 8. In order to 
calculate the NIR number counts, we must fix the UV spec- 
tra. Although we have not included spectral synthesis model 
of stellar populations (but we have used the stellar model 



by Schaerer 2002 in § 2.5.1), the following constant spectral 
energy distribution can be used as a first approximation to 
the NIR number counts: 



L^(M, 2) = Luv(Mvir = M, 



2)/A!AjV : 



(40) 



where Avw is the typical width of the frequency range (~ 
lO^'^-Sx lO^'^ Hz; 400-3000 A in wavelength) where OB stars 
dominate the luminosity. Here we assume Ai^uv = 7 x 10^^ 
Hz. We use the data listed in Table |^ for Luv as a function 

of (Mvir, 2vir). 

In Fig. M, we show the number counts in NIR. The 
solid and dotted lines show the counts for 5 < 2 < 20 and 
5 < 2 < 7, respectively. The vertical dot-dashed line shows 
a proposed detection limit of NGST at the NIR (~ 0.1 njy; 
Gardner & Satyapal 2000). We see that about 2 x 10^ galax- 
ies should be detected per square degree. The detection limit 
can be larger (1 nJy; Mather & Stockman 2000). In this case, 
about 10^ galaxies are detected per square degree. These 
numbers are much larger than those derived for ALMA. In 
other words, all the ALMA sample for 2 > 5 is detected by 
NGST. Therefore, NGST can be used for the two aims: for 
an identification of ALMA sources and for an exploration 
of the universe deeper than ALMA. NGST can also detect 
galaxies whose redshift is larger than 7. The high-redshift 
2 > 5 sample brighter than 0.1 nJy (1 nJy) contains 14% 
(1.2%) of "extremely" high-redshift 2 > 7 sources. 



4.2.3 bSFH vs. cSFH in sub-mm number counts 

Here, we consider the bSFH. The difference between bSFH 
and cSFH is larger in sub-mm number counts than in NIR 
number counts, because dust accumulation makes the FIR 
luminosity increase efficiently while UV suffers from extinc- 
tion by dust. The FIR luminosity increases by an order of 
magnitude owing to the increase of star formation rate and 
dust content. Therefore, the effect of the bSFH can be exam- 
ined by adopting a FIR luminosity 1/10 times smaller than 
for cSFH for 2 > 2burst- For 2 > 2burst, we adopt the same 
star formation rate as that of cSFH. As a result, the FIR 
luminosity rises by 10 times at 2 = 2burst in the bSFH. The 
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Figure 6. Number counts for high-redshift galaxies in the near 
infrared. The solid and dotted lines show the number counts in 
the redshift range of 5 < z < 20 and 5 < z < 7, respectively. The 
vertical dotted line indicates an expected detection limit (0.1 nJy) 
of the iVGST photometry (Gardner & Satyapal 2000). 



expected sub-mm number counts for the bSFH is shown by 
the dashed lines in Figure |^. A clear difference between the 
two scenarios cannot be seen at the ALMA detection limit. 



Since ALMA is not sensitive to galaxies at 2 < 7 (§ 4.2.1), 
we cannot distinguish the two SFHs as long as the difference 
appears at 2 < 7. 

4.3 Integrated light of UV and FIR emissions 

We also predict the level of the fiux integrated for all the 
galaxies from 2max = 20 to 2inin ~ 5, to see whether 
our model is consistent with the current observational con- 
straints, and to compare our result with future more sen- 
sitive observations. The integrated light from all the ex- 
tragalactic objects (extragalactic background radiation) has 
been detected in a wide range of wavelength (e.g., Hauser 
& Dwek 2001). We can examine what fraction of the extra- 
galactic background light is produced by the high-redshift 
galaxies by using our results. 

Theoretically, the intensity of integrated light produced 
by sources between 2inin and 2niax at an observed frequency 
f is estimated by 



/•Zmax /-Afmax 


dz 


dM 




(1 + 2)L(i+2)^ 




47Tdl \ 





dn{M, 2) dV{z) 



dM 



dz 



(41) 



The same spectra (L^) as § 4.2 are adopted for both FIR 
and UV. 

The sub-mm and NIR integrated intensities of the high- 
redshift galaxies are shown in Fig. |^ (thick solid lines). Both 
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Figure 7. Predicted intensity levels of integrated light from all 
the high-redshift (5 < 2; < 20) galaxies for (sub-)millimetre and 
near-infrared wavelengths (thick solid lines). The contribution 
from galaxies with 5 < 2: < 7 is also calculated (thick dotted 
lines). We have also calculated the integrated light from all the 
high-redshift galaxies in the bSFH scenario, but the results are 
indistinguishable from the thick dotted lines in both UV and (sub- 
)millimetre. The observational data for extragalactic background 
radiation are also shown because our predictions should lie be- 
low them. The references are as follows: Bernstein et al. (2002; 
0.8 /im); Cambresy et al. (2001; 1.25 and 2.2 /im); Wright (2001; 
1.25 and 2.2 /xm); Gorjian et al. (2000; 2.2 and 3.5 /xm); Wright 
& Reese (2000; 2.2 and 3.5 ^tm); Hauser et al. (1998; 140 and 240 
/im); Lagache et al. (2000; 140 and 240 /im). We also show the 
ranges estimated from integration of NIR galaxy counts by Totani 
et al. (2001; thick error bars without symbols). The shaded area 
indicates the observed level constrained by COBE measurement 
(extragalactic background; Fixsen et al. 1998). 



fluxes are well below the observed extragalactic background 
radiation. The largest (~ 10%) contribution of the high- 
redshift galaxies is expected in the millimetre range. This 
indicates that future ALMA 1.3 mm observations will be 
the most efficient to isolate the high-redshift contribution. 

In Figure ^, we also show the integrated intensities in 
the case of Zmax ~ 7 and Zmin = 5 (thick dotted lines) . Com- 
paring thick solid and dotted lines, we see that roughly 52% 
and 35% of the high-redshift (2 > 5) integrated hght comes 
from sources with z > 7 at millimetre and UV wavelengths, 
respectively. The model predictions in the case of the bSFH 
are also shown but the results are indistinguishable from 
the thick dotted lines (the difference is at most 7%). There- 
fore, if the bSFH is correct, the integrated hght from 2 > 5 
sources is reduced by half in the millimetre range. 

For NIR, a more elaborate model focusing on Pop 
III stars has been developed by Santos, Bromm, & 
Kamionkowski (2002) and Salvaterra & Ferrara (2002). 
These authors took into account the formation of Pop III 
stars in halos less massive than our Mmin. They also consid- 
ered stars whose mass is larger than 100 Mq because Pop 
III stars are widely believed to be massive. Our model, on 
the other hand, has included Mmin to exclude galaxies that 
cannot sustain gas against energy input from SNe II, but 
adopted a more standard IMF. The integrated hght of NIR 



by Santos et al. (2002) is two orders of magnitude higher 
than our integrated light of high-redshift galaxies, although 
we should note that "extreme" conditions (star formation 
efficiency of 40%, activation of star formation in halos less 
massive than Mmin) are required in their paper to realise 
such a high value of integrated flux. 

4.4 Integrated metal-line intensity 

Suginohara et al. (1999) have proposed that the spatial fluc- 
tuations of integrated metal-line intensity can be used as a 
tracer of metal production in high-redshift sources. Now that 
we have analysed the metal production and the star forma- 
tion history in a consistent manner, it is worth reexamining 
the integrated metal-line intensity within our framework. 

We consider an observation which is sensitive to the fre- 
quency range [i/obs-A!/obs/2, !/obs+Ai/obs/2]. Ifwe are to de- 
tect a metal line whose frequency is i^unc, the redshift range 
which we can observe is [21, Z2] = [i'iine/(t'obs + A:/obs/2) — 
1, fiinc/(^'obs — A:/obs/2) — 1]. Therefore, the integrated line 
intensity, /iino(i'obs; Ai/obs), is 



/linc(i-'oba; AUohs) 



PZ2 r^'In 




/ / 


dM 






dV{z) Lline 




dz 4:Tvd1 


1 


Av is much wider 



dn{M, 



dM 



(42) 



width of each galaxy. In this paper, we assume zi = 5 and 
02 = 7 to conce ntrate on high-redshift galaxies detected by 
ALMA (§ 4.2.1). Here, we again approximate (M, z) with 



(Mvir, Zvir). In Table |, we show I^^^^ relative to the CMB 
intensity, where the superscript "max" indicates that the in- 
tegrated line intensity is calculated by using the maximum 
line intensity as listed in Table ^. We also present Uohslunc^ 
in Table ^ 



As discussed in § 2.5.2, the metal-line intensity should 



be multiplied by < 1 if the gas density of galax ies is lower 



2.1.1 



we obtain 

4 



than n^. Adopting the mean density in 
.7^ ~ 5.5 X 10~^ 9.7 X 10~^ 5.8 x 10'^ 2.9 x 10"", and 
9.7 X 10~^ at 2 ~ 6 (average between zi and Z2) for C I 609 
pm, C I 370 /im, O l 63.2 ^m, O l 146 /xm, and C 11 158 /im, 
respectively. Hence intensities are reduced by the small 
values and are difficult to detect. However, if we adopt, for 
example, the typical density for photo-dissociation region 
in the nearby universe ~ 10^ cm~^ (HoUenbach & McKee 
1979), we obtain ~ 0.68, 0.26, 2.1 x 10'^, 1.0 x 10"^, and 
0.26, for C I 609 ^m, C I 370 fim, O I 63.2 mn, O I 146 
pm, and C 11 158 fim, respectively. In Figure Bl, we show 
as a function of hydrogen number density. Observations 
of the integrated line intensity will constrain the typical gas 
density in high-redshift galaxies by coinparison with theo- 
retical maximum values listed in Table W. If more than two 
metal lines are detected, a consistency check of gas density 
is possible. 



5 SUMMARY AND DISCUSSIONS 

5.1 Summary of evolutionary properties 

In order to quantify the importance of dust on the first star 
formation activity in the universe, we have solved the time 
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Table 4. Metal lines considered in this paper 
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rmax / t 
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(cm 3) 




(nW m~ ) 


C 


0.17 


'Pi- 
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7.93 X 10" 


8 


4.7 X 10^ 


3.6 X lO"-"^ 


4.3 X 10"* 






'P2- 




370 


2.68 X 10" 


7 


2.8 X 10^ 


1.5 X lO"** 


4.6 X 10-5 


O 


1.2 


'Pl- 




63.2 


8.95 X 10- 


5 


4.7 X 10^ 


11 


13 






'Po- 
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1.70 X 10" 


5 


9.5 X 10* 


7.2 X 10"^ 


0.33 


C+ 


0.17 


^Pl/2 - 


- ^P3/2 


158 


2.36 X 10- 


6 


2.8 X 10^ 


1.0 X lO-'' 


5.3 X 10-3 



evolution of dust mass in the galaxies formed in the redshift 
range z > 5, when the age of the universe is < 1 Gyr. We 
have taken into account the importance of H2 abundance 
for the star formation rate, and the formation of molecules 
on dust in a consistent manner (§ ^). In particular, we have 
made the first attempt to tie the star formation efficien cy t o 
H2 abundance in relatively primordial environments (§ 2.2). 
Even when this inefficient phase of star formation is in- 
cluded, an active phase of star formation takes place after a 
few tcii (much shorter than the Hubble timescale) because a 
significant amount of dust is accumulated to activate the H2 
formation on the grain surfaces (Fig. ^). This suggests that 
the grains play an essential role in causing the first active 
phase of star formation. As a result, we have provided a ro- 
bust support for some theoretical works that have implicitly 
assumed that stars are formed at high redshift as efficiently 
as at low redshift. 

Radiative properties of high-redshift star-forming 
galaxies are also predicted. We have found that a significant 
amount of luminosity is radiated in the FIR range. The FIR 
luminosity becomes comparable to the UV luminosity in a 
few tcir, when a significant amount of dust is accumulated 
(Fig. This efficient reprocessing of the UV light into the 
FIR results partly from the dense (i.e., large optical depth) 
environment of high-redshift galaxies. Observations of the 
FIR light (sub-mm light in the observer's restframe) as well 
as those of the stellar light are thus crucial to trace the whole 
stellar radiative energy from high-redshift galaxies. 

In the framework of our model we have also given an 
approximate estimate of some sub-mm metal-line luminosi- 
ties. However, this can only been seen as an upper limit 
to the actual luminosity beca use co nsiderable uncertainty is 
present on the gas density (§ 2.5.2 and Appendix H). If fu- 



ture sub-mm or millimetre observations detect metal lines, 
a density probe of ISM of the high-redshift galaxies will be 
possible. If more than two types of metal-lines are detected, 
density can be estimated more precisely. 



5.2 Future observational tests 

In about ten years, it will become possible to detect sources 
at high redshift {z > 5) in both sub-mm and NIR. The 
luminosity level of the galaxies at these wavelengths will 
put important constraints on our model. Therefore, we have 
calculated the number counts for both wavelengths. As a 
result, we have found that ALMA (450 fim, 850 ^m, and 

1.3 mm bands) and NGST (NIR bands) can detect several 
times 10^ and 10® high-redshift galaxies per square degree, 
respectively. These numbers can be used to test our model. 

We have also calculated the integrated intensity of the 



metal-line emission from the galaxies from z = 5 to 7. Al- 
though a precise determination of the line intensity requires 
a model for the gas density (perhaps a model for photo- 
dissociation region as HoUenbach & McKee 1979 is also nec- 
essary), we can estimate a maximum intensity for the inte- 
grated metal-line intensity. The results are listed in the last 
two columns of Table ^. The contamination with the cosmic 
sub-mm and microwave background could represent a po- 
tential problem, but such high-redshift galaxies have a cor- 
relation scale of the order of 10" (Appendix ^). Therefore, 
a fluctuation analysis of sky brightness in (sub-)millimetre 
can discriminate the metal-line signal from other contam- 
inating sources by examining the typical correlation scale. 
A quantitative analysis of the fluctuations using the struc- 
ture formation theory is left for future work. It may be also 
possible to probe the gas density of high-redshift galaxies 
through the factor T in Appendix ^, if more than two kinds 
of metal-lines are detected. The relative intensities of two 
lines can be used to derive a probable value of tih. 

Finally, we should mention how to select high-redshift 
galaxies efficiently. Galaxy colors (flux ratios between two 
bands) are often used for the selection. We have shown that 
ALMA can detect galaxies at z < 7. The peak of the dust 
emission lies roughly between 500 /im and 800 /im for galax- 
ies between 2 = 5 and 7. For example, 450 /.im vs. 1.3 mm 
flux ratio (450-1300 /im color) gives us a useful informa- 
tion on the redshift, because the peak of flux lies between 
these bands only for the high redshift galaxies. In Figure ^ 
we show the flux at 850 /xm (Sgso Mm) and 450-1300 /im 
color (5450 /im/SLs mm) predicted by the modified black- 
body spectra (eq. p9[ ). The galaxies detected by ALMA in 
that redshift range has typical flux levels ~ 10-100 /iJy and 
the 450-1300 ^.m color is typically less than 3. In Figure ^ 
we also show the ALMA detection limit (horizontal dashed 
line). Galaxies with Mvir > lO" '^ Mq will be detected by 
ALMA. At z > 7, however, the number of such massive 
galaxies is negligible and does not contribute to number 
counts. Typically, galaxies whose redshift is less than 4 fall to 
the right of the vertical dotted line (Takeuchi et al. 2001b). 

We can also select high-redshift galaxies efficiently from 
optical observations by using the "dropout" technique (Stei- 
del et al. 1996). The Lyman limit at the wavelength of 912 
A in the restframe of a galaxy is redshifted to 5500-7300 A 
for galaxies at z = 5-7. Therefore, optical/NIR observations 
of galaxies by NGST (Mather & Stockman 2000) or other 
sensitive facilities provides us a way to sample the high- 
redshift candidates independent from the ALMA sample. A 
large sample of galaxies with V- or J?-band dropout should 
be collected by future observations. After spatial cross iden- 
tification of drop-out sample with ALMA sample, we can 



16 H. Hirashita & A. Ferrara 




1 10 

1.3 mm 

Figure 8. Relation between the flux at A = 850 fira (Ssso fim) 
and the 450-1300 /^m color (S450 ij,m/S\.3 mm)- (We adopt the 
same notations as Fig. 10 of Takeuchi et al. 2001b in this figure.) 
Our model predictions are identified by filled squares for 2 = 5, 
6, 7, and 8 and Mvir = 10", 10^^ ^nd lO^^ Mq. The horizontal 
dashed line indicates the detection limit of ALMA. Low-redshift 
galaxies (2 < 4) typically fall on the right of the vertical dotted 
line (Takeuchi et al. 2001b). 
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Figure 9. Luminosities as a function of formation redshift 2vir 
for Mvir = 10^ '^■^ AfQ. Such a massive galaxy will be detected by 
ALMA. The solid, dotted, and dashed lines represent ultravio- 
let, far-infrared, and O I 146 ^m (maximum; L^'^q) luminosities 
defined at 4 circular times. 



5.3 Connection to lower redshift 

As shown in Figure ^, our predictions connect smoothly to 
the lower-redshift star formation history. Our model, how- 
ever, cannot be applied to galaxy evolution at z < 5 because 
after that epoch dust is supplied from late-type stars as well 
as SNe II. It is observationally known that mergers between 
giant galaxies significantly contribute to luminous infrared 
populations at the local universe (Sanders & Mirabel 1996) 
and even at z ~ 1 (Roche & Eales 1999). When we apply 
our framework to lower redshifts, therefore, it is necessary 
to extend our model to include the details of the merging 
history of galaxies. We should note that the enhancement of 
molecular formation is also a key to star formation activity 
in mergers (e.g., Walter et al. 2002). 

Recent studies using the Subaru telescope (Ouchi et al. 
2002) have pushed observations as deep as z; ~ 5. Therefore, 
the luminosity function (or comoving star formation rate) 
derived from the "Subaru Deep Field", which is as wide as 
600 arcmin^ and as deep as 26 AB magnitude around 7000 
Awill allow us to directly compare our results at z ~ 5. The 
luminosity function at 2: ~ 5 is also important to constrain 
the evolutionary scenario of Lyman break populations found 
at 2 ~ 3. Is the luminosity function of galaxies at z ~ 5 
explained by the same population of Lyman break galaxies 
at 2 ~ 3? Recently, Ferguson, Dickinson, & Papovich (2002) 
have given a negative answer to this question, but further 
studies are necessary to reveal the link between these two 
epochs. 

We have stressed the importance of dust on the for- 
mation of molecular-rich environment. In the lower-redshift 
{z < 5) universe, it is observationally known that there is a 
correlation between the abundances of dust and molecules 
for DLAs (Ge et al. 2001). This strongly suggests the impor- 
tant role of dust for molecule formation (see also Levshakov 
et al. 2002). However, the correlation is not firmly assessed 
and further observational sample seems to be required (Pe- 
titjean et al. 2000). Petitjean et al. also noted that most 
of DLAs may arise selectively in warm and diffuse neutral 
gas. Liszt (2002) has shown that even in a cool medium H2 
formation can be suppressed because of low dust content 
and strong UV radiation field. Moreover, DLA trace a di- 
verse population with various mass, surface brightness, etc. 
(e.g., Pettini 2002). In spite of those complexities, DLAs are 
promising objects to study the link between the abundances 
of dust and molecules in the early universe. 
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APPENDIX A: Ha FORMATION ON GRAINS 

The production rate of molecular fraction via dust surface 
reaction is estimated as (HoUenbach & McKee 1979) 



dt 



= /ofidustTra vS , 



(Al) 



where /hj is the molecular fraction of hydrogen (equation 
|l|), fo is the neutral fraction of hydrogen, nn is the number 
density of hydrogen nuclei, v is the mean thermal speed of 
hydrogen, ndust is the number density of dust grains, a is 
the radius of a grain (spherical grains with a single radius 
are assumed), and S is the sticking efficiency of hydrogen 
atoms. The thermal speed is given by (Spitzer 1978) 



8 ksT 
n mil 



1.4 X 10 



VlOO KJ 



1/2 



(A2) 



where ks is the Boltzmann constant, T is the gas tempera- 
ture, and rriH is the mass of a hydrogen atom. Here, we define 
the reaction rate of the H2 formation on grains, -Rdust, as 

SmuvS 



Rdus 



8a5 
1.4 X 10" 

S 



T 



100 K 



1/2 



0.03 fim 



2 g cm" 



3 -1 
cm s 



(A3) 



where 5 is the mass density of a grain. The dust-to-gas mass 
ratio D can be estimated as 

47ra^(5 Wdust 



3nHmH 



(A4) 



Using -Rdust and V, the right-hand side of equation (Al) 
becomes the second term of equation ([lo|). We take a = 
0.03 /im (Todini & Ferrara 2001). The sticking coefficient 
S is uncertain. Simply, we adopt S ~ 0.2, the value for 
the gas temperature when star formation occurs (T < 300 
K) (HoUenbach & McKee 1979). When T > 300 K, S is 
assumed to be zero. The sticking efficiency of HoUenbach & 
McKee (1979) indicates that the dependence of S on dust 
temperature is negligible as long as the dust temperature, 
Tdust, is less than ~ 75 K. 



from u to I and that from I to u, the fraction of the popula- 
tion in the upper level, /„, is estimated to be 



fu 



Aui + {'yui + 'yiu)n 



(Bl) 



If the coUisional excitation (or deexcitation) occurs on 
shorter timescale than the spontaneous emission, the frac- 
tion in the upper level is 



J U 



lul + llu 



(B2) 



Because of the spontaneous emission, the fraction in the 
upper level is reduced by the factor of T = fu/ fu - For the 
excitation coefficients, the following relation holds (Spitzer 
1978): 



(B3) 



Jlu 

-Jul gi \ ksTc, 

where Qu and gi are statistical weights for the lower and 
the upper levels, respectively, Eui is the energy gap between 
the two states, fce is the Boltzmann constant, and Tex is 
the excitation temperature of the colliding species. In this 
paper, we consider metal lines whose yiu/^ui is order unity 
(i.e., Eui < ksTcx) and whose exciting species is hydrogen 
(n ~ nn). For the metal lines of interest, we finally find that 



nn 



(B4) 



gives a good approximation (correct within a factor of ~ 2), 
where Tier is the critical density defined as 



nil 



Jul 



(B5) 



The line intensity is expected to be ~ L'^nc (§ 2.5.2) if nn is 
higher than (i.e., JT ~ 1). On the contrary, if the density 
if much smaller than n^, the line intensity is quite reduced 
(i.e., T<. 1). 

In Fi gur e Bl, we show as a function of nn, where 



equation (B4) is applied. If we detect line emission intensity 
of high-redshift galaxies, we can make a rough estimate of 
the gas density by comparing the observed line intensity 
with LJfn^. The line ratio between two types of lines can 
also constrain the gas density. 



APPENDIX B: SIMPLE PRESCRIPTION TO 
INCLUDE THE LOW-DENSITY EFFECT ON 
LINE INTENSITY 

Here we describe a simple approximate treatment to calcu- 
late FIR-sub-mm metal-line intensities. Our formula given 
here is appropriate for a one-zone treatment as our model. 
A more accurate treatment requires detailed modeling of 
photo-dissociation regions (HoUenbach & McKee 1979; Tie- 
lens & HoUenbach 1985). 

We consider a population with two energy levels (the 
upper and the lower are labelled as u and I, respectively). 
We consider spontaneous transition, whose Einstein coeffi- 
cient is denoted as A^i, and coUisional excitation by species 
whose number density is n (the coUisional excitation rate 
is expressed as 7!„n) and coUisional deexcitation by the 
same species (the coUisional deexcitation rate is expressed 
as 'yuin). Assuming the equilibrium between the transition 



APPENDIX C: TYPICAL ANGULAR SIZE OF 
CLUSTERING 

The integrated light from galaxies should trace their spatial 
fluctuation produced by gravitational clustering. Therefore, 
if we measure the typical correlation angular scale of the in- 
tegrated light, we can test if the integrated light really comes 
from galaxies. Here, we calculate the typical correlation an- 
gular size of galaxies at a certain redshift. 

The typical mass within a clustered region, Md, is esti- 
mated to be the mass scale corresponding to a{Mci)D{z) = 
1, where D{z) is the growth factor of perturbations (-D(O) = 
1), and a{M) is the variance of density field smoothed with 
mass scale M. In this relation, Md corresponds to the mass 
scale at the break of the Press-Schechter mass function. We 
adopt a{M) given in Appendix of Kitayama & Suto (1996). 
Once Aid is obtained, we can estimate the typical diameter 
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Figure Bl. Reduction factor to the maximum luminosity 
^Unc nictal lines as a function of the number density of liydro- 
gen atoms, nj]- The solid, dotted, dashed, and dot-dashed lines 
represents of C l 609 ^tm, C l 370 /im, O l 63.2 fim, and O I 
146 //m, respectively. Since the critical density of C ii 158 jum is 
the same as that of C I 370 /nm, .7^ of C II 158 jum is represented 
by the dotted line. If the number density nn is known, we can 
obtain the line intensity as J^L'^^^. Inversely, if we obtain by 
comparing observed line luminosity and theoretical L^^, we can 
estimate nn- 



of gravitationally clustered region, Dci, as 

1/3 



3Mci 



(CI) 



where pco is the critical density of the universe at z = (i.e., 
pcof^M is the mean mass density of the universe at z = 0). 

In this paper, we are particularly interested in z ~ 6. 
Using L»(6) = 0.21, we obtain Ma = 2 x 10^° Mq for the 
cosmology assumed in this paper. This corresponds to the 
typical comoving size of 0.33 Mpc (i.e., Dd = 0.05 Mpc). 
This corresponds to an angular scale of 10". 



REFERENCES 

Abel, T., Anninos, P., Zhang, Y., Norman, M. L. 1997, NewA, 2, 
181 

Abel, T., Bryan, G. L., Norman, M. L. 2002, Science, 295, 93 
Anders, E., & Grcvcsse, N. 1989, Gcochim. Cosmochim. Acta, 53, 
197 

Armus, L., Matthews, K., Neugebauer, G., Soifer, B. T. 1998, 

ApJ, 506, L89 
Barkana, R. 2002, NewA, 7, 85 

Bernstein, R. A., Preedman, W. L., Madore, B. F. 2002, ApJ, 571, 
56 

Bromm, V., Coppi, R S., & Larson, R. B. 2002, ApJ, 564, 23 

Buat, v., Dcharvcng, J. M., & Donas, J. 1989, A&A, 223, 42 
Cambrcsi, L., Reach, W. T., Bcichman, C. A., & Jarrctt, T. H. 

2001, ApJ, 555. 563 
Carroll, S. M., Press, W. H., & Turner, E. L. 1992, ARA&A, 30, 

499 

Ciardi, B., Ferrara, A., Governato, F., & Jenkins, A. 2000, MN- 

RAS, 314, 611 
Ciardi, B., & Loeb, A. 2000, ApJ, 540, 687 



Cole, S., Laccy, C. G., Baugh, C. M., & Frcnk, C. S. 2000, MN- 
RAS, 319, 168 

Cox, A. N. 2000, Allen's Astrophysical Quantities (4th ed.; New 
York: Springer) 

Dale, D. A., Helou, G., Neugebauer, G., Soifer, B. T., Prayer, D. 

T., & Condon, J. J. 2001, AJ, 122, 1736 
Devriendt, J. E. G., & Guiderdoni, B. 2000, A&A, 363, 851 
Draine, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 
Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 
Dwck, E., ct al. 1983, ApJ, 274, 168 
Dwck, E., ct al. 1998, ApJ, 508, 106 
Edmunds, M. G. 2001, MNRAS, 328, 223 

Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H., Franceschini, 
A., Fadda, D., & Chary, R. R. 2002, A&A, 384, 848 

Ferguson, H. C, Dickinson, M., & Papovich, C. 2002, ApJ, 569, 
L65 

Ferrara, A. 1998, ApJ, 499, L17 

Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 319, 

539 

Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, C. L., & Shafer, 

R. A. 1998, ApJ, 508, 123 
Frenklach, M., & Feigelson, E. 1997, in ASP Conf. Ser. 122, From 

Stardust to Planetesimals, ed. Y. J. Pendleton & A. G. G. M. 

Tielens (ASP: San Francisco), 107 
Gain, D., & Palla, F. 1998, A&A, 335, 403 
Gardner, J. P., & Satyapal, S. 2000, AJ, 119, 2589 
Gc, J., Bcchtold, J., & Kulkarni, P. 2001, ApJ, 547, LI 
Gehrz, R. D. 1989, in lAU Symp. 135, Interstellar Dust, ed. L. 

J. AUamandola & A. G. G. M. Tielens (Dordrecht: Kluwer), 

445 

Gispert, R., Lagache, G., & Puget, J. L. 2000, A&A, 360, 1 
Gorjian, V., Wright, E. L., & Chary, R. R. 2000, ApJ, 536, 550 
Granato, G. L., Lacey, C. G., Silva, L., Bressan, A., Baugh, C. 

M., Cole, S., & Frenk, C. S. 2000, ApJ, 542, 710 
Haiman, Z., Rees, M. J., & Loeb, A. 1996, ApJ, 467, 522 
Hauscr, M. G., ct al. 1998, ApJ, 508, 25 
Hauscr, M. G., & Dwek, E. 2001, ARA&A, 39, 249 
Hirashita, H., Hunt, L. K., & Ferrara, A. 2002a, MNRAS, 330, 

L19 

Hirashita, H., Tajiri, Y. Y., & Kamaya, H. 2002b, A&A, 388, 439 
HoUenbach, D. J., & McKee, C. F. 1979, ApJS, 41, 555 
Hu, E. M., Cowie, L. L., McMahon, R. J., Capak, P., Iwamuro, 
F., Kncib, J.-R, Maihara, T., & Motohara, K. 2002, ApJL, 

568, L75 

Hutchings, R. M., Santoro, F., Thomas, P. A., & Couchman, H. 
M. P. 2002, MNRAS, 330, 927 

Ikeuchi, S. 1988, Fundam. Cosmic Phys., 12, 255 

Inoue, A. K., Hirashita, H., & Kamaya, H. 2000, AJ, 120, 2415 

Jones, A. P., Tielens, A. G. G. M., & HoUenbach, D. J. 1996, 
ApJ, 469, 740 

Kamaya, H., & Silk, J. 2002, MNRAS, 332, 251 

Katz, N. Furman, I., Biham, O., Pirronello, V., & Vidali, G. 1999, 
ApJ, 522, 305 

Kauffmann, G., & Chariot, S. 1998, MNRAS, 294, 705 

Kcnnicutt, R. C, Jr. 1998, ApJ, 498, 541 

Kitayama, T., & Ikeuchi, S. 2000, ApJ, 529, 615 

Kitayama, T., Susa, H., Umemura, M., & Ikeuchi, S. 2001, MN- 
RAS, 326, 1353 

Kitayama, T., & Suto, Y. 1996, ApJ, 469, 480 

Kozasa, T., Hasegawa, W., & Nomoto, K. 1991, A&A, 249, 474 

Lagache, G., Haffner, L. M., Reynolds, R. J., & Tufte, S. L. 2000, 
A&A, 354, 247 

Lcvshakov, S. A., Dcssaugcs-Zavadsky, M., D'Odorico, S., & Mo- 

laro, P. 2002, ApJ, 565, 696 
Lisenfcld, U., & Ferrara, A. 1998, ApJ, 496, 145 
Liszt, H. 2002, A&A, 389, 393 

Madau, P., Ferguson, H. C, Dickinson, M., Giavalisco, M., Stei- 
del, C. C, & Fruchter, A. 1996, MNRAS, 283, 1388 



Dust and early galaxy evolution 19 



Madau, P., Pozzetti, L., Dickinson, M. 1998, ApJ, 498, 106 
Madau, P., Ferrara, A., & Rees, M. J. 2001, ApJ, 555, 92 
Mather, J. C, & Stockman, H. S. 2000, The Institute of Space 
and Astronautical Science Report SP No. 14, Mid- and Far- 
Infrared Astronomy and Future Space Missions, ed. T. Mat- 
sumoto & H. Shibai (Sagamihara: ISAS), 203 
Malkan, M. A., & Stecker, F. W. 2001, ApJ, 555, 641 
Matsuda, T., Sato, H., & Takeda, H. 1969, Prog. Theor. Phys., 
42, 219 

McKee, C. F. 1989, in lAU Symp. 135, Interstellar Dust, ed. L. 
J. AUamandola & A. G. G. M. Tielens (Dordrecht: Kluwer), 
431 

McKee, C. F., & Ostriker, J. P. 1977, ApJ, 218, 148 
Mo, H. J., fc \Y hitc. S. D. M. 2002, MNRAS, submitted ( [astro 
ph/( f202393| ) 

Moseley, S. H., Dwek, E., Glaccum, W., Graham, J. R., Loewen- 

stein, R. F., & Silverberg, R. F. 1989, Nature, 340, 697 
Nakamura, F., & Umemura, M. 2002, ApJ, 569, 549 
Nagashima, M., Totani, T., Gouda, N., & Yoshii, Y. 2001, ApJ, 
557, 505 

Nishi, R., & Susa, H. 1999, ApJ, 523, L103 

Norman, C. A., & Spaans, M. 1996, ApJ, 480, 145 

Oh, S. P. 2002, MNRAS, submitted 

Omukai, K. 2000, ApJ, 534, 809 

Omukai, K., & Nishi, R. 1998, ApJ, 508, 141 

Ouchi, M., et al. 2002, in RevMexAA Conf. Ser., Galaxy Evolu- 
tion: Theory and Observations, in press 
Partridge, R. B., & Peebles, P. J. E. 1967, ApJ, 148, 377 
Pearson, C. P. 2001, MNRAS, 325, 1511 

Peebles, P. J. E. 1980, The Large Scale Structure of the Universe 

(Peinceton; Princeton University Press) 
Pei, Y. C., Fall, S. M., & Hauser, M. G. 1999, ApJ, 522, 604 
Petitjean, P., Srianand, R., & Ledoux, C. 2000, A&A, 364, L26 
Pettini, M. 2001, in 17th lAP Astrophysics Colloquium, Gaseous 
Matter in Galaxies and Intergalactic Spece, ed. R. Felet et al. 
(Frontier Group: Paris), 315 
Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 
Puget, J.-L., Abergel, A., Bernard, J. -P., Boulanger, F., Burton, 

W. B., Desert, F.-X., & Hartmann, D. 1996, A&A, 308, L5 
Rana, N. C., & Wilkinson, D. A. 1986, MNRAS, 218, 497 
Ripamonti, E., Haardt, F., Ferrara, A., & Colpi, M. 2002, MN- 
RAS, in press 
Roche, N., & Eales, S. A. 1999, MNRAS, 307, 111 
Salvaterra, R. , & Ferrara, A. 2002, in preparation 
Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 



Santos, M. R., Bromm. V.. fc Ka mionkowski. M. 2002, MNRAS, 

submitted (astro-ph/0111467) 
Schaerer, D. 2002, A&A, 382, 28 

Small, I., Ivison, R. J., Blain, A. W., & Kneib, J.-P. 1998, ApJ, 
507, L21 

Shanks, T., Metcalfe, N., Fong, D., McCracken, H., Campos, A., 
& Thompson, D. 2001, in lAU Sump. 204, The Extragalactic 
Infrared Background and its Cosmological Implications, ed. 
M. Harwit & M. G. Hauser (ASP: San Francisco), 347 

Shapiro, P., & Kang, H., ApJ, 318, 32 

Soifer, B. T., Neugebauer, G., Franx, M., Matthews, K., Illing- 
worth, G. D. 1998, ApJ, 501, L171 

Somerville, S. R., & Primack, J. R. 1999, MNRAS, 310, 1087 

Spitzer, L., Jr. 1978, Physical Processes in the Interstellar 
Medium (New York: Wiley) 

Steidel, C. C., Adelberger, K. L., Giavalisco, M., Dickinson, M., 
& Pettini, M. 1999, ApJ, 519, 1 

Steidel, C. C, GiavaUsco, M., Pettini, M., Dickinson, M., & Adel- 
berger, K. L. 1996, ApJ, 462, L17 

Suginohara, M., Suginohara, T., & Spergel, D. N. 1999, ApJ, 512, 
547 

Takeuchi, T. T., Ishii, T. T., Hirashita, H., Yoshikawa, K., Mat- 

suhara, H., Kawara, K., & Okuda, H. 2001a, PASJ, 53, 37 
Takeuchi, T. T., Kawabe, R., Kohno, K., Nakanishi, K., Ishii, T. 

T., Hirashita, H., & Yoshikawa, K. 2001b, PASP, 113, 586 
Tan, J. C, Silk, J., & Balland, C. 1999, ApJ, 522, 579 
Tegmark, M., Silk, J., Rees, M. J., Blanchard, A., Abel, T., & 

Palla, F. 1997, ApJ, 474, 1 
Tielens, A. G. G. M., & HoUenbach, D. 1985, ApJ, 291, 722 
Tinsley, B. M. 1980, Fundam. Cosmic Phys., 5, 287 
Tinsley, B. M., & Danly, L. 1980, ApJ, 242, 435 
Todini, P., & Ferrara, A. 2001, MNRAS, 325, 726 
Tosi, M., & Diaz, A. I. 1990, MNRAS, 246, 616 
Totani, T., & Takeuchi, T. T. 2002, ApJ, 570, 470 
Totani, T., Yoshii, Y., Iwamuro, F., Maihara, T., & Motohara, 

K. 2001, ApJ, 550, L137 
Walter, F., Weiss, A., Martin, C, & Scoville, N. 2002, AJ, 123, 

225 

White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 

Wilson, C. D., Scoville, N., Madden, S. C, & Charmandaris, V. 

2000, ApJ, 542, 120 

Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 

Wright, E. L. 2001, ApJ, 553, 538 

Wright, E. L., & Reese, E. D. 2000, ApJ, 545, 53 

Xu, C, Lonsdale, C. J., Shupe, D. L., O'Linger, J., & Masci, F. 

2001, ApJ, 562, 179 



